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Goals 

The  long-term  goals  are  to  understand  how  the  marginal  ice  edge  zone 
(MIZ)  interacts  with  the  open  ocean  in  an 

oceanographic/meteorological/biological  sense.  This  includes  understanding 
the  processes  and  dynamics  for  the  generation,  dissipation  and  role  of 
mesoscale  eddies  and  chimneys  in  the  MIZ  (and  in  the  open  ocean  after  the 
ice  has  retreated  or  melted)  in  the  Greenland  Sea/Fram  Strait.  A  further 
goal  is  to  understand  how  these  physical  features  interact  with  the 
biological  and  chemical  oceanography  causing  enhancement  of  biological 
production  in  the  MIZ  and  adjacent  open  ocean. 

All  of  this  is  aimed  at  understanding  the  physical  and  biological 
controls  on  the  visible  and  acoustic  environment  of  the  MIZ.  It  is  also 
aimed  at  understanding  the  effects  of  interannual  variability  and  global 
change  which  result  in  variability  in  flux  of  carbon,  through  primary 
production,  to  the  deep  ocean  possibly  buffering  CO^  in  the  atmosphere. 

Objectives 

Our  first  objective  was  to  construct  numerical  models  of  air-ice-ocean 
interactions.  Our  initial  model  (Niebauer  and  Smith,  1989)  is  working  well. 
It  is  being  used  to  help  us  understand  and  to  test  hypotheses  concerning  the 
physical  and  biological  oceanographic  interactions  that  have  been  observed 
as  the  ocean  and  atmosphere  interact  with  the  edge  of  the  Arctic  sea  ice 
pack  in  the  Fram  St ra it -Greenland  Sea.  It  is  also  being  used  to  model  the 
Bering  Sea  MIZ.  The  model  includes  the  interaction  of  ocean  currents  and 
eddys  with  the  ice  edge  as  well  as  the  effects  of  winds,  storms  and  solar 
radiation . 

Our  hypothesis  is  that  the  physical  processes  are  a  key  to  both  the 
biological  phenomena  of  spring  blooms  and  to  bloom  variability. 
Understanding  these  physical  and  biological  phenomena  are  critical  in 
understanding  the  visible  and  acoustic  regime  of  the  marginal  ice  edge  zone. 
Our  objective  is  to  be  able  to  model  these  processes  and  their  variability 
and  ultimately  to  predict  the  phenomena  of  bio-physical  interaction  at  the 
edge  of  the  polar  ice  caps. 

A  second,  fully  3-D,  time  dependent  model  of  the  Greenland  Sea  is  under 
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construction  (Maslowski,  1992)  .  The  objective  here  is  to  model  the  3-D 
aspects  of  the  MIZ  (including  the  biology  and  ice)  but  to  extend  the  space 
scales  to  basin-wide  as  well  as  extending  the  simulations  to  interannual  and 
global  change  time  scales  using  super  computer  resources.  We  have  made 
progress  as  outlined  below,  Weislaw  Maslowski  is  doing  his  PhD  building 
this  model. 

The  second  objective  of  this  contract  was  to  gather  and  analyze  data 
from  the  marginal  ice  edge  zone  to  test  these  hypotheses.  The  data  were 
gathered  from  the  M/V  Polarbjorn  (April-May,  1989)  during  the  Coordinated 
Eastern  Arctic  Experiment  (CEAREX) .  The  scientific  crew  for  this  part  of 
the  pro ject/contract  included  H.J.  Niebauer,  Thomas  Weingartner  and  Steve 
O)c)conen .  The  analyses  of  these  data  continues  with  sections  of  the  data 
published  as  outlined  below. 

Methods 

We  employ  two  basic  methods  to  accomplish  out  objectives;  1)  Numerical 
modeling,  and  2)  A  field  program  from  the  M/V  Polarbjorn  collecting 
coincident  physical  and  biological  oceanographic  data  during  spring  1989 
CEAREX  program  in  the  Fram  Strait  marginal  ice  edge  zone. 

The  numerical  modeling  was  started  in  1985  and  the  basic  model  was 
completed  and  published  in  1989.  The  detailed  methods  are  outlined  in 
Niebauer  and  Smith  (1989) .  The  model  has  been  continuously  improved  upon 
with  support  from  this  contract. 

Basically  the  model  is  made  up  of  three  interacting  sub-models:  1) 
physical  oceanography  (including  meteorology),  2)  biological  oceanography 
and  3)  ice.  The  physical  submodel  is  a  series  of  time-  and  space  dependent 
coupled  non-linear  equations  of  motion  which  include  vertical  and  horizontal 
advection  and  diffusion,  density  gradients,  mesoscale  eddys  as  well  as  wind 
and  insolation  forcing. 

The  ice  submodel  provides  ice  movement,  convergence  and  divergence,  and 
melting  and  freezing. 

The  biological  submodel  includes  two  separate  nutrient  pools  (nitrate 
and  ammonium)  ,  nutrient  upta)ce  by  two  separate  phytoplanlcton  size  classes, 
passive  sin)cing  of  phytoplankton,  grazing  by  both  micro-  and 
macrozooplankton  on  each  phytoplankton  size  class,  export  via  fecal-pellet 
flux  and  material  transfer  to  higher  trophic  levels,  and  nutrient  cycling 
and  regeneration  by  microzooplankton,  macroplankton  and  bacterioplankton . 

The  3-D  physical  model  (with  detailed  bottom  topography)  is  running  in 
the  barotropic  mode  driven  by  historical  winds.  Work  is  progressing  on  the 
baroclinic  mode  using  historical  hydrographic  and  meteorological  data.  Work 
is  starting  on  the  ice  model.  The  biological  modeling  will  begin  within  the 
year  but  will  lean  heavily  on  previous  (Niebauer  and  Smith,  1989)  work. 

Our  methods  for  the  Polarbjorn  field  data  have  been  to  combine  our 
analyses  of  the  physical  (temperature,  salinity  and  density  vs.  depth), 
biological  (chlorophyll  or  phytoplankton)  and  chemical  (nitrate,  ammonium 
etc.)  data  gathered  on  the  CEAREX  biological  cruise  with  the  modeling  with 
the  aim  of  understanding,  first,  the  physics  of  the  eddys,  stratification. 
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etc.,  and  second,  the  physical  controls  or  effects  on  the  biological 
production . 

Results 

Numerical  Modeling 

We  have  constructed  a  numerical  model  of  the  physical  and  biological 
oceanographic  interaction  at  marginal  ice  edge  zones  as  outlined  above.  Our 
experiments  have  emphasized  the  Fram  Strait  in  the  East  Greenland  Sea  in 
spring  (May  and  June)  when  ice  is  actively  melting  and  solar  radiation  is 
increasing,  leading  to  spring  blooms.  However,  the  model  is  generic  and  we 
have  also  applied  it  to  the  Bering  Sea  MIZ  (Niebauer,  1991) .  While  much  of 
the  following  has  resulted  from  previous  ONR  support,  some  has  been  a 
result  of  the  present  support  extending  the  experiments  in  both  the 
Greenland  Sea  and  Bering  Sea. 

We  have  published  several  sets  of  experiments  on  the  Fram  Strait  MIZ 
(Niebauer  and  Smith,  1989;  Smith  and  Niebauer,  In  Press)  including  1)  a 
standard  case  with  oceanic  stratification  and  ice  melt  stratification  but  no 
winds  or  eddys,  2)  variations  of  the  standard  case  with  various 
stratification  schemes  with  and  without  ice  cover,  3)  cases  with  wind 
forcing  and  ice  cover  but  no  eddys,  and  finally,  4)  a  case  with  a  cyclonic 
eddy  under  the  ice  cover. 

For  the  first  two  sets  of  experiments,  whenever  there  is  stratification 
there  is  an  accelerated  accumulation  of  phytoplan)cton  biomass  and  associated 
reduction  in  nitrate.  When  the  stratification  is  due  to  ice  melt  only, 
the  bloom  is  trapped  between  the  ice  and  the  ice-edge  frontal  structure  and 
stratification  as  is  observed  in  field  studies.  When  there  is  additional 
stratification  away  from  the  ice  edge,  the  phytoplanlcton  biomass  is  spread 
out  but  remains  spatially  correlated  with  the  stratification.  Regardless  of 
the  source  of  stratification,  the  magnitude  and  timing  of  the  bloom  is 
similar.  That  is,  maximum  chlorophyll  coincides  with  minimum  nitrate  (this 
occurs  in  10  days  to  two  wee)cs)  but  lags  maximum  primary  production  by  two 
to  three  days.  The  stratification,  which  enhances  the  bloom,  also  restricts 
the  introduction  of  nutrients  from  below  and  limits  primary  production  after 
about  one  wee)c .  Nitrate  levels  fall  from  8-10  uM  to  less  than  1  uM  while 
chlorophyll  levels  increase  from  less  than  1  ug  1"'  to  over  5  ug  1 
Zooplan)cton  biomass  increases  in  response  to  the  bloom  reaching  a  maximum 
about  two  wee)t3  after  maximum  phyt op lani^ton  biomass.  Both  the  maximum 
primary  production  rates  (>1  gm  C  m’^  day”')  and  the  maximum  zooplan)cton 
biomass  (>22  mg  dry  wt  m'^)  are  in  quantitative  agreement  with  field 
estimates.  In  the  absence  of  stratification,  a  slow  increase  in 
chlorophyll  occurs  but  no  nutrient  limitation  is  observed.  The  nanoplanlcton 
and  microzooplankton  do  not  play  a  significant  role  in  these  experiments. 

In  the  third  set  of  experiments,  with  wind  forcing,  the  major  finding 
was  that  regardless  of  which  way  the  wind  stress  was  directed,  in  either  the 
Fram  Strait  or  the  Bering  Sea,  there  was  enhancement  of  primary  production. 
When  the  wind  stress  was  directed  on-ice  or  parallel  to  the  ice  edge  with 

the  ice  to  the  right  of  the  wind,  some  ice  edge  upwelling  was  observed  while 

when  the  wind  was  off-ice  or  parallel  to  the  ice  with  the  ice  edge  to  the 

left  of  the  wind,  some  downwelling  was  observed.  However,  the  actual 

upwelling  or  downwelling  advection  had  a  minor  effect  on  the  primary 
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production.  The  reasons  for  this  are  that  the  region  of  vertical  advection 
is  narrow  and  localized  to  the  ice  edge  which  was  moving  in  response  to  the 
winds.  In  addition,  the  relatively  strong  melt  water  derived  density 
stratification  resists  vertical  movement.  For  the  case  of  ice  being  blown 
over  open  water  (the  downwelling  cases) ,  the  most  noticeable  effect  was  that 
the  ice  e.xerts  a  drag  on  the  surface  of  the  water  carrying  a  nitrate-rich 
layer  seaward.  After  the  winds  cease  there  was  a  region  of  about  12  1cm 
seaward  of  the  ice  edge  that  was  rich  in  nitrate  and  melt-derived 
stratification,  and  enhanced  growth  resulted.  However,  productivity  at  the 
initial  ice  edge  decreased  due  to  increased  ice  cover. 

In  the  cases  where  the  ice  was  blown  off  the  ocean  (the  upwelling 
cases)  the  most  noticeable  effect  was  to  increase  the  light  available  for 
photosynthesis  in  recently  ice  covered  regions  that  were  replete  with 
nutrients  leading  to  increased  productivity.  Again,  upwelling  seemed  to 
have  little  effect  on  the  hydrographic  structure  and  thus  had  little  effect 
on  the  productivity  at  the  ice  edge. 

There  is  a  noticeable  difference  between  the  circulation  associated 
with  upwelling  situations  and  circulation  associated  with  downwelling.  The 
difference  is  due  to  the  way  the  ice  reacts  to  wind  stress.  When  the  ice  is 
forced  out  past  the  ice  edge  over  open  water  (the  downwelling  situation), 
the  ice  tends  to  thic)cen  or  "pile  up"  down  wind  and  so  the  effect  of  wind 
stress  on  the  ocean  via  the  ice  is  enhanced.  When  the  ice  is  forced  back 
into  the  icepack  (the  upwelling  situation)  the  ice  tends  to  thin  at  the  ice 
edge.  With  lower  concentrations  of  ice,  less  effect  of  the  wind  is  felt  in 
the  ocean  and  so  the  mesoscale  ice  edge  ocean  circulation  is  not  so  strong. 
In  the  model  the  upwelling  velocities  were  1/3  the  downwelling  velocities. 
Biologically,  the  thicker  ice  covering  more  of  the  ocean  restricted  light  in 
the  photic  zone.  When  the  ice  was  thinned  by  the  wind,  productivity  was 
more  wide  spread  due  to  the  increase  in  light. 

In  the  last  of  the  four  initial  cases,  a  cyclonic  eddy  was  spun  up 
under  the  ice  cover  for  a  week  and  then  allowed  to  relax.  (We  did  not  do 
these  experiments  for  the  Bering  Sea  because  eddys  have  not  yet  been 
observed  up  on  the  shelf  at  the  MIZ.)  During  spin  up  there  was  upwelling  in 
the  center  of  the  eddy  with  surface  divergence  which  caused  a  divergence  in 
the  ice  cover  which  cleared  the  ice  cover  over  the  eddy.  The  strongest 
horizontal  currents  were  at  the  surface  at  >20  cm  s”'  with  flow  away  from 
the  center  of  the  eddy  at  0.25  km  day”'.  Maximum  upwelling  velocities  were 
40  m  day"'  below  50  m  in  a  core  about  6  km  wide.  This  core  is  flanked  by 
columns  of  downwelling.  This  circulation  gave  the  hydrographic  structure 
the  appearance  of  upwelling  in  the  upper  layers  but  downwelling  below  100  m. 
Similar  structure  has  been  seen  in  nitrate  cross  sections  from  the  Fram 
Strait . 

When  the  eddy  forcing  was  turned  off  there  was  downwelling  in  the 
center  of  the  eddy  and  surface  convergence  but  not  enough  to  cause  the  ice 
to  recover  the  surface  of  the  ocean.  We  estimated  that  the  eddy  spin-down 
would  take  about  two  months  (ie.,  the  conversion  of  potential  energy,  stored 
in  the  upwelled  density  structure  during  eddy  spin-up,  back  to  kinetic 
energy  of  circulation  via  Coriolis  dynamics) . 

The  eddy  circulation  swept  the  ice  off  surface  waters  that  were  high  in 
nitrate  concentrations.  This  caused  a  sudden  increase  in  the  light  regime 
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which,  when  combined  with  the  nutrient  rich  water,  led  to  increased 
productivity.  The  timing  and  magnitudes  of  maximum  primary  productivity, 
nitrate,  chlorophyll  and  zooplankton  were  nearly  identical  with  those  of  the 
first  two  sets  of  experiments.  The  only  real  difference  was  that  the  bloom 
was  trapped  between  two  ice  edges  and  that  the  stratification  under  the  open 
water  was  reinforced  by  upwelling. 

The  experiments  in  the  Bering  Sea  were  basically  quite  similar  to  those 
in  the  Fram  Strait  with  some  exceptions.  Nutrient  levels  in  the  Bering  Sea 
are  twice  those  in  the  Fram  Strait  so  typically  the  primary  production  is 
higher  (10  mg  m~^  vs.  as  high  as  40  mg  m"^  in  the  Bering  Sea) .  While 
nutrient  levels  play  a  role,  this  difference  is  also  due  to  the  differences 
in  environment  between  a  shelf  (Bering  Sea)  and  a  deep,  open  ocean  (Fram 
Strait)  system.  The  depth  of  water  at  the  Bering  Sea  MIZ  is  much  more 
shallow  than  the  Fram  Strait  (100  m  vs.  3000  m)  so  that  the  Bering  Sea  water 
is  colder  (less  oceanic)  during  the  spring  MIZ  bloom  which  tends  to  exclude 
zooplankton  (too  cold)  so  that  chlorophyll  levels  are  higher  (no  grazing) . 
Much  of  this  production  goes  to  the  shallow  benthic  communities.  In  the 
Fram  Strait  MIZ,  because  the  water  is  more  oceanic,  the  average  temperature 
tends  to  be  warmer  in  the  water  column  so  there  are  more  zooplankton  biomass 
grazing  the  phytoplankton  so  there  is  less  chlorophyll  in  the  water.  So  in 
the  Fram  Strait  we  are  beginning  to  look  less  at  passive  sinking  of 
phytoplankton  to  take  carbon  from  the  marine  atmosphere  to  the  deep  ocean, 
and  beginning  to  look  more  at  modeling  aggregation  of  phytoplankton 
processes  to  ship  the  primary  production  past  the  zooplankton  to  the  deep 
ocean,  as  well  as  the  more  well  known  fecal  pellet  sinking. 

Field  Experiments 

A  main  objective  of  our  contract  was  to  gather  and  analyze  the  ADCP 
data  from  the  M/V  Polarbjorn  BIO  cruise  in  the  Fram  Strait.  We  used  the  ONR 
support,  as  well  as  NSF  support,  to  develop  ADCP  processing  and  analysis 
methods.  We  applied  this  analysis  initially  to  a  set  of  ADCP  data  from  the 
MIZ  in  the  Bering  Sea.  Horizontal  currents  were  acquired  with  an  ADCP  along 
5  cross-sections  through  the  Bering  Sea  MIZ  in  April-May,  198'7.  Tidal 
current  calculations,  using  a  coincident  current  meter  mooring  for 
calibration,  were  subtracted  from  the  ADCP  data  to  extract  the  residual 
horizontal  currents.  Vertical  currents  were  then  calculated  through 
continuity  in  which  along-ice  variability  was  considered  to  be  negligible. 
Vector  plots  of  the  vertical  circulation  reveal  circulation  cells  within 
w)iich  the  most  energetic  upwelling  occurs  on  the  unstratified  side  of  the 
frontal  structure  and  the  most  energetic  downwelling  occurs  on  the 
stratified  side  of  the  frontal  structure. 

In  support  of  the  observed  circulation,  the  maximum  in  the  mean 
integrated  chlorophy 1 1 / nit rogen  ratio  from  these  same  5  cross-sections 
coincides  with  upwelling  while  the  minimum  in  the  mean  chlorophyll/nitrogen 
ratio  coincides  with  downwelling.  The  positive  correlation  between  the  mean 
biochemical  activity  and  mean  vertical  flow  corroborates  the  existence  of 
persistent  and  relatively  well  defined  vertical  circulation  cells  associated 
with  the  Bering  Sea  MIZ  frontal  structure  as  well  as  suggesting  that  the 
technique  of  using  ADCP  to  calculate  vertical  flow  is  valid. 

The  great  disappointment  of  our  work  was  the  failure  of  the  ADCP  on  the 
Polarbjorn  during  the  CEAREX  BIO  cruise.  We  saw  many  features  in  the  Fram 
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strait  MIZ  from  the  hydrographic  data  as  outlined  below.  In  light  of  our 
work  with  the  ADCP  data  from  the  Bering  Sea,  I  believe  that  we  would  have 
made  significant  progress  had  the  ADCP  not  broken  down. 

There  are  additional  interesting  features/processes  that  show  up  in  the 
CEAREX  BIO  cruise  data  that  are  being  analyzed,  modeled  and  written  up  (eg.. 
Smith  et  al.,  1991).  For  example,  chlorophyll  concentrations  were  found 
quite  deep  in  the  water  column  (~100-200  m) ,  which  is  below  the  photic  zone, 
that  were  still  taking  up  nitrate  or  still  photosynthetically  viable.  We 
originally  thought  that  this  must  be  the  effect  of  sinking  in  chimneys  but 
closer  examination  of  the  physical  data  suggests  that  this  may  be  the  result 
of  downwelling  in  anti-cyclonic  eddies  (or  alternatively,  downwelling 
between  2  cyclonic  eddies).  The  data  do  not  show  a  clear  cut  chimney  with 
sinking  in  the  chimney.  The  data  do  show  strong  sinking/downwelling  in  what 
look  like  anticyclonic  features.  The  energetics  of  these  eddies  (eddy 
streets)  must  still  be  sorted  out  and,  in  fact,  progress  is  being  made  on 
the  eddy  and  eddy  streets  as  results  from  the  3-D  model  outlined  above. 
Particle  aggregation  is  also  being  investigated. 

We  also  get  the  impression  from  the  data  that  chimneys  may  form  as  a 
result  of  cyclonic  eddy  circulation  that  forces  salty  but  somewhat  warmer 
water  to  the  surface  where  atmospheric  interaction  can  draw  out  more  heat. 
The  result  is  cold,  salty  water  at  the  surface  which  then  sinks. 
Preliminary  results  show  that  the  model  can  handle  this  type  of  process. 
Our  aim  here  is  to  continue  to  do  model  runs  to  test  ideas  on  whether  the 
chimneys  and/or  the  eddys  ate  causing  downwelling  of  primary  production,  or 
some  combination.  An  additional  complication  is  that  if  ice  is  present,  ice 
may  insulate  the  ocean  from  cooling  and  convection  and  further,  if  the  ice 
is  melting,  the  surface  water  may  too  fresh  to  sink. 

Work  continues. 
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